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Out-of-Band  Response  of  Reflector  Antennas 


David  A,  Hill 

Electromagnetic  Fields  Division 
National  Bureau  of  Standards 
Boulder,  CO  80303 


The  response  of  reflector  antennas  to  out-of-band  frequencies  has 
been  analyzed  using  physical  optics.  A simple  approximate  ex- 
pression has  been  obtained  for  the  effective  aperture,  and  this 
expression  yields  both  the  receiving  pattern  and  the  frequency 
dependence  of  the  on-axis  gain.  The  theory  has  been  compared  with 
published  out-of-band  measurements,  and  the  pattern  agreement  is 
good,  but  the  measured  gain  falls  below  the  theory.  This  discrep- 
ancy is  caused  by  mismatch  loss  in  the  coax-to-waveguide  adapter, 
and  this  mismatch  loss  has  been  analyzed  theoretical ly.  The  basic 
physical  optics  model  for  symmetrical  reflectors  has  been  extended 
to  include  offset  and  dual  reflectors,  reflector  surface  rough- 
ness, and  transient  excitation. 

Key  words:  coax-to-waveguide  adapter;  effective  aperture;  focal 
region;  paraboloid;  physical  optics;  Poynting  vector;  out-of-band 
response;  receiving  pattern;  reflector  antenna. 


1.  Introduction 


The  response  of  antennas  to  out-of-band  frequencies  [1-3]  plays  an  impor 
tant  role  in  interference  and  jamming  problems.  Reflector  antennas  are  of 
particular  interest  because  they  are  used  so  frequently  and  because  they  have 
a strong  response  to  above-band  frequencies.  The  analysis  of  reflector  anten 
nas  at  above-band  frequencies  is  complicated  by  the  presence  of  higher-order 
modes  which  can  propagate  in  typical  waveguide  feeds.  Frequencies  well  below 
the  in-band  frequency  are  not  important  because  they  are  cut  off  by  typical 
waveguide  feeds.  Consequently,  "out-of-band"  will  refer  only  to  above-band 
frequencies  throughout  this  report. 

An  earlier  out-of-band  analysis  of  reflector  antennas  [4]  yielded  the 
radiation  pattern  for  each  propagating  mode  and  performed  a statistical  analy 
sis  based  on  the  statistics  of  the  higher-order  mode  coefficients.  The  probl 
em  is  that  the  statistics  of  the  higher-order  modes  are  highly  dependent  on 
the  specific  feed  system  and  are  not  generally  known.  Such  problems  can  be 
avoided  by  analyzing  reflector  antennas  in  the  receiving  mode  and  computing 
only  the  total  received  power  carried  by  all  the  propagating  modes.  This 
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total  power  is  easily  related  to  a generalized  effective  aperture  which  ap- 
pears to  be  the  most  convenient  receiving  characteristic  for  out-of-band  fre- 
quencies. Both  the  receiving  pattern  and  the  frequency  response  are  given  by 
the  effective  aperture,  and  the  results  depend  only  on  the  antenna  param- 
eters. The  power  which  is  coupled  from  the  waveguide  to  the  detector  does  of 
course  depend  on  the  details  of  the  feed  system,  but  that  portion  of  the  prob- 
lem can  be  analyzed  separately.  In  any  case,  the  total  power  in  the  waveguide 
is  a useful  upper  bound  for  the  detected  power. 

The  organization  of  this  report  is  as  follows'.  Section  2 contains  an 
analysis  of  the  fields  and  the  Poynting  vector  in  the  focal  region  of  a sym- 
metrical paraboloid  for  plane-wave  incidence.  In  Section  3,  the  Poynting 
vector  is  integrated  over  the  aperture  of  the  feed  horn  to  yield  the  total 
received  power.  The  integration  can  be  done  numerically  or  analytically,  and 
a simple  analytical  approximation  is  adequate  for  most  cases.  The  theory  is 
compared  with  some  earlier  out-of-band  measurements  for  frequencies  from  3 GHz 
(in  band)  to  10  GHz  [5],  and  the  agreement  is  fairly  good.  Section  4 contains 
a number  of  extensions  to  the  theoretical  model,  such  as  offset  and  dual  re- 
flector geometries,  the  effect  of  surface  roughness,  and  the  response  of  a 
typical  coax-to-waveguide  adapter.  Transient  excitation  is  also  considered. 
Section  5 summarizes  the  results  of  this  study  and  makes  recommendations  for 
future  work . 

2.  Fields  in  the  Focal  Region  of  a Paraboloid 

In  this  section  we  derive  expressions  for  the  electric  and  magnetic 
fields  and  the  Poynting  vector  in  the  focal  region  of  a symmetrical  parabolic 
reflector.  Much  of  the  previous  work  on  the  focal  region  fields  of  para- 
boloids [6-10]  was  directed  toward  design  of  feed  systems,  and  only  on-axis 
incidence  was  considered.  More  recently,  Valentino  and  Toulios  [11]  computed 
the  electric  field  in  the  focal  region  of  an  offset  paraboloid  for  off-axis 
incidence.  For  simplicity  we  consider  only  the  symmetrical  paraboloid,  but  we 
extend  the  previous  derivation  [11]  to  include  the  magnetic  field  and  the 
Poynting  vector  for  off-axis  incidence. 


2 


2.1  Physical  Optics  Integration 

The  perfectly  conducting,  parabolic  reflector  is  shown  in  figure  1.  The 
reflector  has  a diameter  D and  a focal  length  f.  The  origin  of  a rectangular 
coordinate  system  (x,y,z)  is  located  at  the  focus  of  the  paraboloid.  Our 
derivation  and  notation  follow  Valentino  and  Toulios  [11]  fairly  closely,  but 
we  introduce  several  simplifications.  Because  we  consider  the  symmetrical 
paraboloid,  the  offset  angle  0Q  in  [11]  is  zero.  For  the  symmetrical  para- 
boloid, without  loss  of  generality  we  assume  the  incident  plane  wave  is  inci- 
dent in  the  xz  plane  ( 4>s  = 0 in  [11]).  Finally,  we  assume  that  the  field 
point  is  located  in  the  focal  plane  (z2  = 0 or  02  = tt/2  in  [11]). 

The  physical  optics  surface  current  J on  the  reflector  is  given  by 


J = 2 n x H,  , 


(1) 


where  n is  the  unit  normal  to  the  reflector  and  H_.  is  the  incident  magnetic 
field.  For  plane  wave  incidence,  FL  is  given  by 


H. 

-i 


-H 


H.  e 

l 


(2) 


where  k is  the  free  space  wavenumber  ( =2n/X) , X is  the  wavelength,  and  the 
exp  (jwt)  time  dependence  is  suppressed.  In  all  cases,  y with  a subscript 
denotes  a unit  vector.  We  chose  the  quantity  to  be  zero  at  the  center  of 
the  reflector,  and  it  is  given  by 


8 = u • r 

s -p  -s 


(3) 


where  rc  is  directed  from  the  center  of  the  reflector. 


An  infinitesimal  surface  current  patch  of  area  dS  at  a point  P on  the 
reflector  produces  electric  and  magnetic  fields  [8,11]  at  the  point  ?2: 

-jkR 


- Ur  (J  • Ur)]  • 

ik  e"jkR 

‘ 17  u * Hr)  — K—  dS  - 


(4) 
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where  n is  the  intrinsic  impedance  of  free  space.  As  indicated  in  figure  1,  R 
is  the  distance  between  P1  and  P2,  and  uR  is  directed  from  P to  P2.  Also  it 
is  assumed  that  kR  is  much  greater  than  unity. 


In  order  to  evaluate  eq  (4),  we  restrict  our  analysis  to  the  region  near 
the  focus  (r2«r1).  Thus  we  can  assume  that  uR  - -u^  and  R - r1  except  in 
the  phase  term.  The  required  expressions  for  the  unit  normal  n and  the  sur- 
face differential  dS  have  been  given  by  Bern  [12]: 


and 


3-ri 

n = x 


30 


3 r 3 r 3 r 

-1  , -1  -1 

/ |—  x 


3 <J> 


30 


3 (j)  • 


(5) 


8 El 1 *  3-rl 
dS  = x d8i  d*i  • 


3 (j) 


Substituting  eq  (5)  into  eq  (4)  and  carrying  out  some  of  the  differentiations, 
we  obtain 

dE  = -j  — - e J'k^s+R^(r1  sin01  d<j>1  d0x)  . 
rLx  1 


{[(u  • uM)  + 


1 


3r 


(u  • u..)]  u 


-<)>]_  -H  r1sin01  3^  -r1  -H  -0L 


(6) 


1 

"[(V  U-H>  + 7^<V  U-H)]  V 


and 


dH  = — — e Jk^s+  R^(ij  sin01  dcj>1  d0x) 


rLXn 
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Vtt 
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^ (u-r/  “h)]  -e/  (d, 


* UH)  } » 
l “H  -<p  ^ 


where  Ej  = nH-j.  The  expression  for  dE  agrees  with  previous  results  [11,12]. 
Note  that  dE  and  dH  are  orthogonal  and  both  dE  and  dH  are  transverse  to  u^ 
This  is  consistent  with  the  large  kR  assumption.  The  expression  for  R can  be 
approximated  as  follows 
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(7) 


R = r:  [1  - 2 — sine1  sine2  cos  - 4>2) 


r 2 '2  1/2 

-2  — cos01  cos02  + 2 ] 


1 


= rl  - r2  [si n 0 L sin02  cos(^1-cj)2)  + cos01  cos©2]  . 

In  eq  (7)  we  have  neglected  terms  in  r22/r i,  and  the  validity  of  this  approx- 
imation has  been  discussed  previously  [11].  The  phase  term  ft$  in  eq  (3)  can 
be  written 


ft  = (f  - r , ) cos0  + r,  sin©,  sin0  co s<t> i 
s 1 s 1 1 s 1 


(8) 


Also,  rx  is  given  by 

2 f 

Pl  1 + COS01 


(9) 


The  polarization  of  the  incident  magnetic  field  yH  can  be  written 


u,,  = a u + a u + au 
-H  x -x  y-y  z-z 


where  a2  + a2  + a2  = 1. 
x y z 


(10) 


Without  loss  of  generality,  we  confine  the  analysis  to  two  orthogonal  polari- 
zations. For  the  incident  magnetic  field  polarized  in  the  xz  plane,  we  have 


a = cos0  , a = 0,  a = -sin0  . (11) 

x s y z s 

For  the  incident  magnetic  field  polarized  perpendicular  to  the  xz  plane,  we 
have 


ax  = 0,  ay  = 1,  az  = 0 . (12) 

In  order  to  evaluate  the  integration  on  ^ from  0 to  2n  in  eq  (6),  we 
take  advantage  of  the  fact  that  the  integrands  are  periodic  in  4>1 . The  eval- 
uation of  the  4> L integrations  is  given  in  Appendix  A,  and  the  resultant  ex- 
pressions for  E and  H are 
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(13) 
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J e o 
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xn 
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-z 


zn  zn 

The  upper  limit  em  of  the  0!  integration  is  given  by  [11] 


_i 

0m  = tan  (D/4  f)  . 
m 


(14) 


All  other  quantities  in  eq  (13)  are  defined  in  Appendix  A. 

In  general,  the  0!  integration  must  be  performed  numerically,  and  a com- 
puter code  has  been  written  to  evaluate  eq  (13).  Note  that  the  summations 
truncate  at  n = 2 which  is  in  contrast  to  the  case  of  the  offset  parabola 
[11,12]  where  the  summations  run  to  n = °°  . 
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2.2  Large  f/D  Approximation 


When  f/D  is  large,  then  <3m  as  given  by  eq  (14)  is  small,  and  small  argu- 
ment approximations  can  be  used  for  the  tr igonometric  functions  of  ©i  . In 
this  case,  eq  (13)  simplifies  to 


- Fo  ^ay  -x  ax  -y^  *0m  * 
H “ ~~  (ax  ux  - ay  uy)  I6m  , 


(15) 


where 


E = -2 jkf  E.e’J'kf  sin2  (0  /2) 


5 


I 


0m 


g 

cosec2  (0  /2)  m 

j j ( k 0 x P ) s i n 0 x d0x  , 

2 o 

o 


and  P = / (x2  - f sin0^)?  + y22  . 


By  setting  sin0,  - 0,  and  cosec2  (0  12)  - (2/0  )2 

m m 

ate  I Qm  analytical ly: 

2 (k©  P) 

I 

0m  k © P 

m 


in  eq  (15),  we  can  evalu- 


(16) 


Note  that  E and  H are  orthogonal,  and  both  are  transverse  to  y . The  results 
in  eqs  (15)  and  (16)  are  consistent  with  earlier  approximations  [8,11]. 

A fairly  simple  physical  picture  can  be  obtained  from  eqs  (15)  and 
(16).  The  maximum  of  IQ  occurs  for  P = 0,  and  this  occurs  for  the  following 
focal  plane  coordinates: 


x2  = f sin0^  and  y2  = 0.  (17) 

The  point  determined  by  eq  (17)  is  essentially  the  geometrical  optics  point 

for  a ray  incident  on  the  center  of  the  reflector.  For  0s  equal  to  zero  the 

maximum  fields  occur  at  the  focus,  but  for  0 not  equal  to  zero  the  maximum  is 

s 
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shifted  as  indicated  by  eq  (17).  Away  from  the  maximum,  the  decay  is  more 
rapid  for  higher  frequencies  (larger  k)  and  also  for  larger  reflectors 
(larger  0 ) . Also  the  peak  electric  field  EQ  is  larger  for  higher  frequencies 
and  for  larger  reflectors. 

2.3  Poynting  Vector 

The  Poynting  vector  is  of  particular  interest  because  in  Section  3 we 
will  integrate  it  over  the  aperture  of  the  feed  horn  in  order  to  determine  the 
total  received  power.  The  real  Poynting  vector  S is  given  by 

S = i Re  (E  X H*)  , (18) 

where  Re  denotes  the  real  part  and  * denotes  complex  conjugate.  A computer 
code  has  been  written  to  evaluate  eq  (18)  using  the  integral  expressions  for  E 
and  H in  eq  (13) . 

For  large  f/D,  we  can  substitute  eqs  (15)  and  (16)  into  (18)  to  obtain 
the  following  expression 


2 n y 


+ a.2) 


2J, (ke  P) 
1 m 

k0  P 
m 


u 

-z 


(19) 


For  the  incident  magnetic  field  polarized  in  the  y direction,  a = 1 and  ax  = 

0.  For  the  magnetic  field  polarized  in  the  xz  plane,  a = cosQ  and  a =0. 

x s y 

Since  we  are  interested  in  small  scan  angles,  we  can  replace  cos  by  unity 


and  rewrite  eq  (19)  as 


S = u S 

- -z  o 


2J, (k0  P) 
1 m 

¥1  P 

m 


where  S = - — 
o 2 


Eo'2 


(20) 


Thus  the  approximate  2 is  independent  of  the  polarization  of  the  incident 
field  and  contains  only  a z component  Sz . 
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Most  reflector  antennas  do  not  have  large  f/D  ratios,  but  the  approxima- 
tion in  eq  (20)  turns  out  to  be  surprisingly  good  even  for  relatively  small 
values  of  f/D.  Figure  2 shows  S2  for  three  different  values  of  f/D  for 
the  on-axis  case  0$  = 0.  The  curves  for  f/D  = 1.0  and  0.4  were  computed  from 
the  general  integral  expressions  in  eq  (13),  and  the  curve  for  f/D  = 00  was 
computed  from  eq  (20).  We  see  that  the  dependence  on  f/D  is  quite  weak,  and 
the  large  f/D  approximation  in  eq  (20)  is  adequate  for  most  cases.  For  off- 
axis  incidence  (8  * 0),  eq  (20)  is  also  a good  approximation,  and  figure  3 

shows  some  results  for  f/D  = 1.  The  large  f/D  curves  were  computed  from  eq 
(20),  and  the  f/D  = 1 curves  were  computed  from  (13)  for  f = 1 OX  and  for 
polarization.  The  dependence  on  f/X  and  the  incident  polarization  is  fairly 
weak.  Note  that  the  agreement  with  the  large  f/D  approximation  is  better  for 
the  smaller  scan  angle  0$  = 10°,  but  is  not  too  bad  for  0$  = 20°.  We  are 
not  interested  in  very  large  0$  because  the  entire  physical  optics  method  be- 
comes questionable  for  0$  too  large. 

3.  Received  Power 
3.1  Feed  Horn  Response 

Initially  we  consider  a feed  horn  with  an  aperture  of  arbitrary  shape. 
Normally  feed  horn  aperture  dimensions  are  on  the  order  of  a wavelength  at  the 
in-band  frequency.  Consequently,  the  aperture  dimensions  can  be  assumed  to  be 
electrically  large  at  the  higher  out-of-band  frequencies.  Thus  we  make  the 
simplifying  Kirchhoff  approximation  that  the  fields  in  the  aperture  are  equal 
to  the  incident  fields. 

If  the  incident  fields  are  nonuniform,  as  in  the  focal  region  of  a para- 
boloid, then  the  power  passing  through  the  feed  horn  aperture  Pr  is  given  by 
the  integral  of  the  Poynting  vector  over  the  aperture. 

P = - / S dA  . (21) 

r A 2 

Sz  is  the  z component  of  the  incident  Poynting  vector,  and  the  geometry  is 
shown  in  figure  4.  Since  we  assume  that  no  power  is  dissipated  in  the  walls 
of  the  horn  and  waveguide,  the  power  propagating  down  the  waveguide  is  also 
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given  by  Pr  in  eq  (21).  Because  the  waveguide  will  normally  be  multimoded  at 
out-of-band  frequencies,  Pr  is  the  total  received  power  in  all  the  propagating 
modes.  The  simple  theory  in  eq  (21)  does  not  give  the  individual  waveguide 
mode  amplitudes,  bp,  but  it  gives  the  sum  of  the  square  of  the  modal  ampli- 
tudes, £ | b ^ | 2 , which  is  equal  to  the  total  power  if  the  modes  are  properly 
normal ized. 

A more  detailed  analysis  of  the  feed  horn  and  the  junction  between  the 
feed  horn  and  the  waveguide  would  be  required  in  order  to  determine  the  indi- 
vidual waveguide  mode  coefficients.  We  have  stayed  away  from  this  additional 
complexity  because  a knowledge  of  the  individual  mode  coefficients  is  probably 
not  very  useful.  Any  waveguide  bends,  transitions,  or  irregularities  would 
produce  mode  conversion  and  a change  in  the  mode  coefficients.  Even  in  the 
uniform  waveguide,  the  field  distribution  would  change  along  the  waveguide 
because  the  individual  modes  have  different  phase  velocities.  In  contrast, 
the  total  power  remains  constant  along  the  waveguide. 

For  the  simple  case  where  the  incident  field  is  a uniform  plane  wave, 
then  Sz  is  simply  given  by 


S 


z 


-S  coso  , 
o 


(22) 


where  SQ  is  the  incident  power  density,  and  0 is  the  incidence  angle  shown  in 
figure  4.  In  this  case  the  integral  in  eq  (21)  is  easily  evaluated  to  yield 


Pr  = Pq  cose  , 


(23) 


where  P^  = AS„  . 
o o 

Thus  we  have  the  simple  result  that  the  receiving  pattern  of  an  electrically 
large  receiving  horn  is  simply  cose  and  is  independent  of  polarization  and 
the  detailed  shape  of  the  aperture.  This  result  only  holds  when  we  consider 
the  total  multimode  power  and  even  then  is  a high  frequency  approximation 
which  neglects  edge  diffraction.  If  we  consider  the  receiving  pattern  for  an 
individual  waveguide  mode,  then  the  pattern  has  a lobe  structure  as  expected 
for  an  electrically  large  antenna.  The  specific  case  of  an  open-ended,  paral- 
lel plate  waveguide  is  analyzed  in  Appendix  B in  order  to  illustrate  the  dif- 
ference between  the  mode  patterns  and  the  total  power  pattern. 
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3.2  Circular  Aperture  Integration 


For  a feed  horn  aperture  of  arbitrary  shape,  the  integral  for  the  total 
received  power  in  eq  (21)  must  be  evaluated  numerically.  This  numerical  eval- 
uation can  be  rather  time  consuming,  but  for  the  special  case  of  a circular 
aperture  we  can  obtain  an  analytical  approximation  to  eq  (21). 


For  a circular  aperture  of  radius  pm,  the  integral  in  eq  (21)  can  be 
written 

2 tt  pm 

P r = _ / j S^p2  dp2  d<t>2>  (24) 

o o 

where  p2  = /x2  + yf,  <j>2  = tan  ^(y2/x2), 

and  x2  and  y2  are  the  focal  plane  coordinates  as  shown  in  figure  1.  For  the 
general  case  where  Sz  must  be  determined  by  numerical  integration,  the  p2 
and  <t>2  integrations  must  also  be  done  numerically,  and  a computer  program  has 
been  written  for  this  case.  For  large  f/D,  we  can  substitute  eq  (20)  into  eq 
(24)  and  obtain  the  following: 

2w  pm  2 J i ( k O P ) 2 

Pr  = - S J j [ ke  pm  ] P2  dp2  d*2.  (25) 

oo  m 


First  we  evaluate  eq  (25)  for  axial  incidence  (0$  =0).  In  this  case,  P 
= 92  and  the  <t>2  integration  simply  yields  a factor  of  2 tt: 


m 


P = - 2tt  S 
r o 


2 Ji(kemp2)  2 

[ — \7r~z: — J p 2 dp2. 


kemP2 


(26) 


The  92  integration  can  be  done  to  yield  the  following  result: 


where  P 

o 


P = P ri  - J (ke  p ) 
r o L ov  mm' 


* (D/2)2. 


(27) 


PQ  is  the  total  power  incident  on  the  reflector,  and  the  factor  in  brackets 
varies  from  zero  to  unity  as  the  quantity  ke  p increases  from  zero  to  infin- 
ity.  Both  the  power  density  Sz  and  the  received  power  Pr  are  shown  in  figure 
5.  The  ratio  P r /PQ  is  called  aperture  efficiency,  and  typical  values  for  in- 
band  reception  are  on  the  order  of  50  to  90  percent. 
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The  result  in  eq  (27)  has  been  given  previously  by  Minnett  and  Thomas 
[8],  but  they  did  not  treat  the  case  of  off-axis  incidence  (es*0).  For  the 
case  of  0s*  0,  we  lose  the  4>  symmetry  and  must  resort  to  an  approximate  inte- 
gration method.  Consider  first  the  case  where  f | si n e$ | is  less  than  pm.  Then 
the  maximum  intensity  (or  geometrical  optics  point)  is  located  at  X2  = f sines 
and  y2  = 0 inside  the  circle  p2  = Pm  as  shown  in  figure  6.  In  order  to  do  the 
integration  analytically,  we  change  the  integration  area  from  a circle  of 
radius  pm  centered  at  the  origin  to  two  semicircles  of  radii  p+  and  p_  cen- 
tered at  the  geometrical  optics  point  as  shown  in  figure  6.  In  this  case,  the 
integration  can  be  evaluated  as  follows: 


p+  2 Ji(k0  p)  2 

- ' So  f / l ~Te/-  ) 


p dp 


P-  2 Ji(ke  p) 

I 1 

o m 


(28) 


= — S-  [2  - Jn(k0mpj  - J1(k0mP.)  - J_(k0  P ) - J i (k0  P )], 

2 o m + 1 m +'  o m - 1 m - 

where  p = pm  ± f sin0$.  Note  that  for  0S  = 0,  eq  (28)  reduces  to  eq  (27). 

For  the  case  where  p | s i n 0S | is  greater  than  pm,  we  use  a similar  strategy 
as  indicated  in  figure  7.  In  this  case,  the  circular  area  is  replaced  by  an 
annular  sector,  and  the  received  power  integral  can  be  written 


p+  2 J i ( k e p)  2 

Pn  “ - 2cf>  SA  / m 1 

r Yc  o J L 

p- 

<f> 


“ko  p 
m 


p dp 


(29) 


= P — [J  (k0  p ) + J i ( k 0 p ) - J (k0  PJ  - J i ( k 0 PJ], 
o ir  o'  m - 1 m - o'  m +'  mH+y  ’ 


-1  pm 

where  P+  = f sine$  + p , p.  - f sine$  - pm  and  *c  = sin  ( rTfwr  ). 

s 

It  is  easy  to  see  that  when  f sin0s  is  equal  to  pm,  then  p_  = 0,  4>c  = 
tt/2,  and  eqs  (28)  and  (29)  yield  the  same  result  for  Pp.  Also,  it  can  be  seen 
from  figure  7 that  eq  (29)  provides  an  upper  bound  for  Pr  because  the  sector 
includes  all  of  the  circular  area  and  the  integrand  is  always  positive. 
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3.3  Effective  Aperture 


The  effective  aperture  A0  is  defined  as  the  received  power  divided  by  the 
intensity  of  the  incident  plane  wave  [13] 

Ae  = Pr/(|Ei|2/2n)  = Pp/(-S0).  (30) 


In  our  case  we  use  the  same  definition  but  recognize  that  Pr  is  the  received 
power  contained  in  all  the  propagating  modes.  In  general,  A0  must  be  eval- 
uated numerically  because  Pr  must  be  evaluated  numerically  as  indicated  by  eq 
(21).  However,  for  the  circular  aperture  we  can  use  the  approximate  expres- 
sions for  Pp  given  by  eqs  (28)  and  (29).  In  doing  so,  it  is  convenient  to 
normalize  Ae  to  the  physical  aperture  Ap  = tt  ( D/ 2 ) 2 . Then  we  can  write  the 
normalized  effective  aperture  as 


1 2 2 2 2 

4 [2  - J (ke  PJ  - Ji(ke  PJ  - J (ke  P ) - J,(ke  p )],  f sine„  < P 
2 o'  m +'  iV  m +'  o'  m - iV  m - ’ s n 


= 4 


(31) 


<t>r  2 2 2 2 

— [J  (ke  p ) + J i ( k e p ) - J (ke  p ) - J L ( k e p,)],  f sine„  > p 
tt  o'  m - iV  m - o'  m +'  m +'  s n 


where  p , = If  sine  ± p and  <t>  is  defined  in  eq  (29). 

±i  s mi  Yc  ^ 

It  is  possible  to  study  the  wide  angle  behavior  (f  sines  >>  e ) by  using 
the  large  argument  approximations  for  the  Bessel  functions  [14]  in  eq  (31). 


J (x)  « — cos2x  and  Ji(x)  ~ — sin2x.  (32) 

From  eqs  (31)  and  (32),  we  obtain  the  following  asymptotic  expansion  for  Ae/Ap 

A 4 p 2 

_e m 

A ir2ke  (f  sine  )3’ 
pm  s 

The  (sine$)“3  dependence  in  eq  (33)  is  the  same  as  the  wide  angle  dependence 
for  the  radiation  pattern  of  a circular  aperture  of  constant  illumination 
[13].  The  pm^  dependence  is  to  be  expected  because  it  is  proportional  to  the 
area  of  the  receiving  horn. 
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The  angular  dependence  of  A_/An  for  various  values  of  ke  p is  shown  in 

figure  8.  The  limiting  case  of  k0mPm  = 00  corresponds  to  geometrical  optics 

where  the  power  is  focused  to  a single  spot.  When  (f/pm)  sin0s  is  less  than 

unity,  the  spot  is  inside  the  receiving  horn,  and  Ag/Ap  is  unity.  When  (f/pm) 

sin0s  is  greater  than  unity,  the  spot  is  outside  the  receiving  horn,  and  Ae/Ap 

is  zero.  For  k0  p = 10,  the  focused  spot  is  smeared  out,  and  the  pattern  is 
mm 

widened.  For  k0  p =3,  the  pattern  is  widened  further.  For  (f/pm)  sin0_ 
m m m s 

much  greater  than  unity,  the  patterns  in  figure  8 begin  to  approach  the  asymp- 
totic expansion  given  by  eq  (33). 

An  approximate  expression  for  the  beamwidth  0b  can  be  obtained  by  setting 
0^/2  equal  to  the  incidence  angle  where  the  first  zero  in  the  Bessel  func- 
tion in  eq  (20)  is  located  at  the  edge  of  the  receiving  horn.  Then  the 

Poynting  vector  in  the  receiving  aperture  will  be  small,  and  the  integrated 
power  and  effective  aperture  will  be  small.  This  condition  can  be  written 

f sin(0b/2)  = Pm  + pi,  (34) 

where  JL  ( kp  1 0 ) = 0 and  kp^  = aL  « 3.832  represents  the  first  zero  of  Ji 
[14]. 


The  above  equation  is  easily  solved  for  0^: 


0,  = 2 sin 
b 


1 Pm  + ai/k0m 
-_1  ( -m- -m- ) 


p + a1/k0m 
m 1 m 


(35) 


When  k is  infinite,  we  recover  the  geometrical  optics  result,  sin0^  = 2 Pm/f, 
which  is  shown  in  figure  8.  As  k becomes  smaller,  0b  becomes  larger  as  shown 
in  figure  8.  For  sufficiently  small  k,  we  have 


sin0,  « 


2 <*! 

I <Tf 
m 


4 a] 

To 


2.44X 
D ' 


(36) 


The  above  result  agrees  with  the  beamwidth  between  the  first  nulls  for  a cir- 
cular aperture  with  constant  illumination  [13]. 


The  patterns  in  figure  8 were  determined  from  the  analytical  approxima- 
tion in  eq  (31).  The  accuracy  of  eq  (31)  has  been  compared  with  numerical 
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integration  for  numerous  cases.  A typical  comparison  for  f/D  = 1 is  shown  in 
figure  9.  The  numerical  integration  curve  was  done  for  the  incident  magnetic 
field  polarized  in  the  y direction,  while  the  analytical  approximation  is  in- 
dependent of  the  incident  polarization.  For  the  numerical  integration 
results,  there  is  a slight  dependence  on  polarization.  Note  that  the  general 
agreement  is  fairly  good,  but  that  the  analytical  result  is  somewhat  larger 
for  (f/p  ) sin0s  greater  than  unity.  This  is  primarily  a result  of  the  upper 
bound  feature  of  the  sector  area  approximation  in  figure  7.  The  numerical 
integration  is  time  consuming  because  it  involves  triple  numerical  integra- 
tion, over  Gj  in  eq  (13)  and  over  p£  and  <j>2  in  eq  (24).  Consequently,  it  is 
simpler  and  more  efficient  to  use  the  analytical  expression  in  eq  (31),  and 
this  expression  is  adequate  for  most  cases. 

In  figures  10-13,  we  compare  the  analytical  approximation  with  previously 
reported  experimental  results  [5]  for  frequencies  from  3 GHz  (the  in-band  fre- 
quency) to  10  GHz.  The  antenna  tested  had  a diameter  D of  1.22  m (4  ft)  and 
an  f/D  ratio  of  0.32.  The  small  f/D  ratio  provided  a good  test  for  the  ana- 
lytical approximation  which  is  best  for  large  f/D.  The  feed  horn  was  rectan- 
gular, 7.21  cm  x 5.00  cm  (2.84  in  x 1.97  in)  and  for  the  theoretical  compari- 
son was  modeled  by  a circular  aperture  with  radius  pm  = 3.39  cm  to  yield  the 
same  area.  The  pattern  data  in  [5]  was  taken  over  a range  of  -40  to  +40°  in 
the  H plane  and  over  a smaller  range  (either  -5  to  +5°  or  -10  to  +10°)  in  the 
E plane.  Consequently,  we  show  only  the  H plane  comparisons  in  figures  10- 
13.  The  pattern  comparisons  are  for  relative  power  or  effective  aperture. 
Strictly  speaking,  the  theory  should  not  be  used  at  the  in-band  frequency,  3 
GHz,  in  figure  10  because  the  feed  horn  and  waveguide  are  not  electrically 
large. 

A difficulty  in  comparing  the  out-of-band  measurements  with  the  theory  in 
figures  11-13  is  that  the  coax-to-wavegui de  adapter  has  an  unknown  response 
for  multimode,  out-of-band  excitation.  Our  analysis  in  Section  4.5  indicates 
a very  complicated  response  as  a function  of  frequency  and  modal  content. 
Because  our  theory  yields  the  total  power  in  the  waveguide  and  the  adapter  has 
a small  response  for  some  modes  and  some  frequencies,  we  can  expect  the  theory 
to  provide  an  upper  bound  or  envelope  for  the  detected  power.  This  seems  to 
be  the  case  in  figures  11-13. 
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In  figure  14  we  show  the  measured  gain  for  the  same  antenna  [5],  To  con- 
vert our  effective  aperture  Ae  to  gain,  we  divide  by  \2/4ti,  the  effective 
aperture  for  an  isotropic  antenna  [13],  Also  shown  is  the  gain  result  using 
the  physical  aperture  Ap.  The  effective  aperture  curve  is  seen  to  gradually 
approach  the  physical  aperture  curve  as  the  frequency  is  increased.  In  con- 
trast, the  measured  gain  becomes  highly  variable  as  the  frequency  is  in- 
creased, and  this  is  thought  to  be  due  to  the  irregular  frequency  response  of 
the  adapter.  These  measured  results  point  out  the  importance  of  the  adapter 
in  determining  what  portion  of  the  total  waveguide  power  is  actually  coupled 
into  the  system. 

4.  Extensions  to  the  Model 
4.1  Distant  Sidelobes 

Because  the  results  in  the  previous  section  are  based  on  physical  optics, 
the  patterns  are  not  expected  to  be  valid  for  large  values  of  0$.  Some 
interesting  comparisons  of  measured  patterns  and  physical  optics  calculations 
for  in-band  frequencies  have  been  made  by  Stubenrauch  and  Yaghjian  [15]. 
Typically,  they  found  that  physical  optics  gave  good  accuracy  for  the  main 
lobe,  less  accuracy  for  the  first  few  sidelobes,  and  poor  results  for  the  dis- 
tant sidelobes. 

For  out-of-band  frequencies  we  are  not  able  to  calculate  the  details  of 
the  pattern  accurately  as  indicated  in  figures  11-13,  but  we  find  that  physi- 
cal optics  gives  a reasonable  estimate  of  the  envelope  pattern  over  the  entire 
forward  hemisphere  (0$  < n/2) . To  illustrate  this  point,  we  use  the  geometri- 
cal theory  of  diffraction  (GTD)  to  provide  an  estimate  of  the  edge-diffracted 
field.  In  figure  15,  a plane  wave  is  incident  on  the  top  edge  of  the  reflec- 
tor an  an  angle  0S  from  the  axis,  and  we  are  interested  in  the  diffracted 
field  near  the  focus.  If  the  magnetic  field  incident  on  the  edge  is  , then 
the  diffracted  field  Ha  at  the  focus  is  [16] 

-jkr 

H = H,  ^ £.  D , (37) 

“ 1 / r e 

e-J*/4  6 

where  D = C 

e 2 / 2nk 
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and  C 


1 


1 


cos  1/2  (9  -0.  ) cos  V2T0  . + 9.  )* 

c scat  l nc  c scat  i nc 


In  C,  the  minus  sign  applies  to  the  case  where  the  incident  magnetic  field  is 
polarized  parallel  to  the  plane  of  incidence,  and  the  positive  sign  applies  to 
perpendi cul ar  polarization.  For  completeness,  we  could  add  the  ray  diffracted 
from  the  lower  edge,  but  we  neglect  it  because  we  are  just  obtaining  the  order 
of  magnitude  of  the  diffracted  field.  Also,  we  are  neglecting  the  curvature 
of  the  diffracting  edge.  The  angles  and  distance  required  in  eq  (37)  are 
given  by 


0scat 
and  r 

e 


, e 

. _JH  a 
2 2*  w i nc 


2f 

1 + COS0  * 

m 


(38) 


If  we  assume  that  pm  is  much  less  than  rg,  then  the  power  density  inci- 
dent on  the  aperture  is  constant,  and  the  z component  Sz  is  given  by 


sz  = - \r\  I Hd  |2  cosem. 


(39) 


The  effective  aperture  A0 

A 

e 


is  then  given  by 

- S up2 
z m 

| H | 2 


2 D2  /r  . 
me  e 


(40) 


If  we  normalize  l\e  by  the  physical  aperture  Ap  = u(D/2)^,  and  use  eqs  (37)  and 
(38)  in  (40),  we  obtain 


A p 3 C2  cose  (1  + cos0  )3 
e _ m m m 

A kf3  64u  sin2  0 

p m 


(41) 


It  is  interesting  that  the  asymptotic  physical  optics  expression  (eq  (33))  has 
the  same  factor  p ^/(kf^),  but  a different  angular  dependence.  The  pm^  depen- 
dence is  simply  a result  of  A0  being  proportional  to  the  feed  aperture  area. 
The  k 1 dependence  is  typical  of  diffracted  power,  and  the  f dependence  is  a 
result  of  the  Ap  normalization  (f dependence)  and  the  inverse  distance  (f“M 
dependence  of  the  edge-diffracted  power. 

There  is  a range  of  incidence  angle,  n/2  <0  < n - 0 , where  the  feed 

3 3 * s m’ 

horn  is  also  directly  illuminated  by  the  incident  field.  The  geometry  for 
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this  case  is  shown  in  figure  16.  In  this  case  it  is  easy  to  show  that  the 
normalized  effective  aperture  due  to  direct  illumination  of  the  feed  horn  is 

A 2 p 2 

7T  = ( — [p  ) cos  (»  - e ),  f < e$  < * - e (42) 

P 

In  figure  17,  we  show  the  various  wide  angle  results  for  effective  aper- 
ture for  the  reflector  parameters  given  in  figure  11.  Strictly  speaking,  the 
physical  optics  result  is  not  valid  beyond  em/2  (=  38°)  because  shadowing  of 
the  reflector  takes  place,  but  in  practice  it  is  still  probably  the  best  esti- 
mate out  to  0$  = 90°.  Beyond  90°,  direct  illumination  and  edge  diffraction 
can  yield  a larger  value  of  Ag.  The  direct  illumination  only  extends  to  the 
shadow  boundary  at  0s  = n - 6m  = 104°.  In  the  vicinity  of  the  shadow  bound- 
ary, the  Keller  diffraction  coefficients  as  given  by  eq  (37)  are  not  valid, 
but  can  be  replaced  by  an  integral  form  [16].  However,  the  maximum  diffracted 
field  is  only  one-half  the  incident  field  even  at  the  shadow  boundary.  Note 
that  the  GTD  results  depend  on  the  polarization  of  the  incident  field. 

Although  the  diffraction  coefficients  decrease  as  0$  increases  beyond  120°, 
there  is  the  possibility  of  a backlobe  at  0$  = 180°  because  the  entire  reflec- 
tor rim  can  contribute  an  in-phase  diffracted  field. 

In  summary,  it  is  probably  most  convenient  to  use  the  physical  optics 
expression  for  es  < 90°  and  to  recognize  that  the  effective  aperture  can  in- 
crease somewhat  for  e$  > 90°  because  of  edge  diffraction  or  direct  illumina- 
tion of  the  feed  horn.  To  try  to  obtain  an  accurate  result  for  large  0$  is 
probably  unrealistic  because  the  results  depend  on  the  construction  details 
(struts,  reflector  edge,  etc.)  of  the  antenna.  Based  on  a large  number  of 
radiation  patterns,  the  CCIR  [17]  uses  the  isotropic  value  (A0  = /4tt ) for 
reflector  antennas  at  large  angles  for  interference  calculations. 

4.2  Transient  Fields 

When  a transient  wave  is  incident  on  the  reflector  antenna,  the  waveform 
undergoes  dispersion.  If  the  frequency  spectrum  of  the  incident  waveform  is 
known,  then  the  problem  can  be  analyzed  in  the  frequency  domain,  and  the  final 
result  can  be  transformed  to  the  time  domain.  Another  approach  is  to  obtain 
the  impulse  response  of  the  reflector  and  to  convolve  the  incident  waveform 
with  the  impulse  response.  This  method  has  the  advantage  that  the  impulse 
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response  can  often  yield  a useful  physical  interpretation.  Also,  the  duration 
of  the  impulse  response  provides  an  estimate  of  the  pulse  stretching  which 
occurs  when  the  impulse  response  is  convolved  with  an  incident  pulse.  In  this 
section  we  obtain  a simple  approximate  expression  for  the  impulse  response  of 
the  reflector,  but  do  not  attempt  to  analyze  the  effect  of  the  feed  horn  or 
the  waveguide  feed. 


In  eqs  (15)  and  (16),  the  approximate  expressions  are  given  for  _E_  and  H_ 
in  the  focal  region.  For  the  approximations  made  in  eqs  (15)  and  (16),  the 
pol arizations  of  _E_  and  H_  are  the  same  as  their  incident  polarizations  and  the 
ratio  of  E to  H is  the  free  space  impedance  n.  Thus  it  is  sufficient  to  treat 
the  scalar  problem  for  either  E or  H,  and  the  polarization  factors  in  eq  (15) 
can  be  ignored.  Thus  we  can  write  the  electric  field  spectrum  as 


E(u>) 


2 J i ( k 0 P) 
1 v m ' 

ke  P 
m 


(43) 


where  Eq  = -2 jkf  E^w)  e ^ s i n 2 ( 6^/ 2 ) and  P is  defined  in  eq  (15)  as  the 
distance  from  the  geometrical  optics  point  (x2  = f S1’n8s>  yz  = 0).  E^w)  is 
the  spectrum  of  the  incident  waveform,  and  for  impulse  excitation  E^w)  is 
unity.  By  setting  E^w)  = 1,  making  the  small  0m  approximation,  and  writing 
in  terms  of  the  modified  Bessel  function  Ij_  [3.2],  we  can  rewrite  eq  (43) 


D e 


E(«) 


2P 


1 1( j“b) 


(44) 


where  b = P0  /c. 

m 


The  impulse  response  E(t)  is  the  inverse  Fourier  transform  of  eq  (44) 

E(t)  = F_1[E(u.)]  = \ E ( oi)  eJajtdu.  (45) 


The  inverse  transform  of  I^(jwb)  can  be  found  in  Laplace  transform  (s  = ju>) 
tables  [18] 


F“  [I^jwb)]  = 


- t 


, It  I < b 


it  b / b2  - t2 

0 , It  I > b. 


(46) 
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The  exponential  factor  in  eq  (44)  introduces  a time  delay,  f/c,  because  the 
phase  reference  was  taken  as  the  center  of  the  reflector.  Consequently, 
E(t)  can  be  written 


!(t)  =2^n<T>’ 


(47) 


x/b 

/ (x/b)2  -1 


where  E ( x)  = 
n 


cl  < b 


, |t|  > 0 


and  x = t - f/c. 


Thus  the  approximate  impulse  response  has  a finite  time  width  w given  by 

2 


w = 2b  = 


2 P 0 2 0 / (x9  - f sin0  )2  - y9: 

m m ^ s 


(48) 


The  shape  of  the  normalized  impulse  response  E (x)  is  shown  in  figure  18.  The 
singularities  at  x = +b  are  integrable  square  root  singularities  and  present 
no  difficulty  in  convolving  E(t)  with  any  realistic  incident  waveform. 

When  P is  zero,  both  eqs  (44)  and  (47)  are  indeterminate.  However,  if  we 
take  the  limit  of  P approaching  zero,  eq  (44)  becomes 

EM  “ - TT1  j“  e"J“f/c,  P = 0.  (49) 


The  inverse  transform  of  jw  [19]  is  the  unit  doublet,  6 ' ( t ) , and  the  inverse 
transform  of  eq  (49)  is 


D0 

E(t)  - - 6’(t),  P = 0. 


(50) 


When  the  doublet  eq  (50)  is  convolved  with  an  incident  waveform,  the  incident 
waveform  is  differentiated  [19]. 


From  eq  (48)  we  see  that  the  width  of  the  impulse  response  is  propor- 
tional to  P,  the  distance  from  the  geometrical  optics  point.  For  on-axis 
incidence  (0$  = 0),  the  most  distant  points  of  interest  are  at  the  edge  of  the 
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circular  aperture,  P = pm.  As  an  example,  consider  the  parameters  in  figures 
10-14:  13  = 1.326  and  Pm  = 3.39  cm.  In  this  case,  eq  (48)  yields  a pulse 

width  w = 0.3  ns.  However,  if  we  consider  off-axis  incidence,  then  P can  be 
much  larger.  If  we  again  consider  the  antenna  parameters  from  figures  10-15 
(f  = 39  cm  and  9m  = 1.326),  an  off-axis  incidence  angle  of  0S  = 30°  yields  a 
pulse  width  of  w = 1.73  ns  at  the  center  of  the  aperture  (x2  = y2  = 0). 

It  should  be  stressed  that  the  results  in  this  section  are  based  on 
physical  optics  and  are  not  valid  at  low  frequencies.  This  means  that  the 
transient  results  in  eqs  (47)  and  (50)  are  not  valid  for  large  times. 

However,  if  the  incident  waveform  has  very  little  low  frequency  content,  the 
convolution  results  should  still  be  fairly  accurate.  The  problem  of  disper- 
sion caused  by  the  feed  horn  and  the  waveguide  has  not  been  addressed,  and 
such  results  would  depend  on  the  detailed  geometry  of  the  horn  and  guide. 
Dispersion  in  waveguides  has  been  studied  by  inverse  transforms  of  the  modal 
series  [20]  for  the  acoustic  problem,  and  the  same  technique  would  be  valid 
for  the  electromagnetic  problem. 

4.3  Reflector  Roughness 

The  analysis  in  this  report  has  assumed  a smooth  parabolic  reflector,  and 
this  assumption  is  valid  when  the  tolerance  of  the  reflector  surface  is  suf- 
ficiently small  in  terms  of  the  wavelength.  The  effect  of  surface  roughness 
has  been  studied  in  the  classic  paper  of  Ruze  [21]  and  in  more  recent  theoret- 
ical papers  [22,23].  In  general,  the  problem  is  very  complicated,  but  we  can 
make  use  of  the  simple  expression  from  Ruze  [21].  By  a statistical  analysis 
of  the  physical  optics  integral,  Ruze  derived  the  following  expression  for  the 
radiation  pattern  G(0,<j>): 

G(0,<|>)  - G ( 0 , 4>)  e + G ( 0 , 4>) , (51) 

o s 

where  GQ  is  the  radiation  pattern  of  the  smooth  reflector,  62  is  the  mean 
square  phase  error  over  the  reflector  aperture,  and  Gs  is  a scatter  term  which 
is  normally  negligible  in  the  vicinity  of  the  main  beam.  The  mean  square 
phase  error  is  given  by 


= (4tt£/ A)  2, 


(52) 
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where  e is  defined  as  the  effective  reflector  tolerance.  For  large  f/D,  e is 
approximately  equal  to  the  reflector  tolerance  measured  normal  to  the  surface, 
and  for  smaller  f/D  the  effective  tolerance  e is  somewhat  less  than  the  normal 
surface  tolerance  [21], 


For  cases  where  the  second  term  in  eq  (51)  is  negligible,  the  effect  of 
reflector  roughness  is  simply  given  by  the  exponential  factor  in  eq  (51). 
Thus  we  can  define  a roughness  loss  L (the  ratio  of  GQ  to  G) : 


62  ( 4ne/X )2 

L = e = e 


(53) 


Normally,  L will  be  small  for  in-band  frequencies  because  the  reflector  toler- 
ance e will  be  small  compared  to  the  in-band  wavelength  X.  Figure  19  shows 
curves  of  roughness  loss  as  a function  of  frequency  for  various  values 
of  e/\  where  XQ  is  the  in-band  wavelength.  It  is  seen  that  roughness  loss 
can  increase  rapidly  with  increasing  frequency  because  of  the  exponential  de- 
pendence of  eq  (52).  The  key  parameter  is  the  surface  tolerance  e which  is 
not  normally  known,  but  has  been  inferred  by  Ruze  [21]  for  some  antennas  from 
the  frequency  dependence  of  the  antenna  gain.  Typical  values  of  e are  on  the 
order  of  a millimeter  or  less.  For  e = 1 mm  and  a frequency  of  3 GHz  (xQ  = 

10  cm),  the  ratio  e/\  is  0.01,  and  the  center  curve  in  figure  19  would 
apply.  Thus  the  roughness  loss  would  be  negligible  at  3 GHz,  but  would  be 
approximately  7 dB  at  30  GHz. 

For  large  values  of  e/X,  the  second  term  in  eq  (50)  becomes  dominant,  and 
a very  broad  pattern  results.  This  term  depends  on  a number  of  parameters 
including  the  correlation  distance  of  the  surface  roughness  [21]  and  will  not 
be  considered  here.  However,  it  is  worth  noting  that  this  scatter  term  can 
actually  raise  the  pattern  level  in  the  distant  sidelobes.  Strictly  speaking, 
the  derivations  of  Ruze  [21]  apply  to  the  in-band  case,  but  his  results  should 
also  apply  to  the  effective  aperture  at  out-of-band  frequencies  because  both 
cases  are  described  by  physical  optics  integrals. 

4.4  Offset  Parabolas  and  Dual  Reflectors 

The  analysis  in  Sections  2 and  3 was  done  for  a symmetrical  paraboloid 
with  a prime  focus  feed,  but  other  reflector  antenna  configurations  are  also 
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of  interest.  In  general,  other  configurations  are  more  complicated  to  analyze 
rigorously,  but  in  many  cases  can  be  described  by  an  equivalent  symmetrical 
paraboloid. 

Consider  first  the  offset  paraboloid  shown  in  figure  20.  The  offset 
geometry  has  the  advantage  of  eliminating  blockage  by  the  feed  horn.  The 
focal  region  fields  of  this  geometry  have  been  analyzed  by  Bern  [12]  for  on- 
axis  incidence  (es  = 0)  and  by  Valentino  and  Toulios  [11]  for  arbitrary  inci- 
dence. Typically,  the  focal  region  fields  of  the  offset  paraboloid  are  ap- 
proximately equal  to  those  of  the  symnmetrical  paraboloid  if  the  equivalent 
focal  length  f1  and  equivalent  diameter  D'  are  used  in  place  of  f and  D.  Thus 
the  formulas  in  Sections  2 and  3 are  directly  applicable  to  the  offset  para- 
boloid if  the  equivalent  parameters  are  used.  For  an  offset  angle  0Q,  the 
equivalent  focal  length  is  given  by  [12] 


f' 


2 f 

1 + COS0  * 
0 


(54) 


The  equivalent  diameter  D'  is  simply  the  projected  diameter  in  the  plane  per- 
pendicular to  the  axis  of  the  parabola  as  shown  in  figure  20.  The  equivalent 
parameter  method  is  most  accurate  for  small  values  of  6 but  Bern  [12]  has 
found  good  results  for  0Q  as  large  as  45°.  Also,  Valentino  and  Toulios  [11] 
have  shown  that  the  geometrical  optics  shift  as  predicted  by  eq  (17)  holds  for 
incidence  angles  of  0S  = 5°  and  10°  for  an  offset  parabola  with  0$  = 44°. 


Symmetrical  dual  reflector  antennas  can  also  be  analyzed  by  replacing 
them  with  single  reflectors  of  equivalent  focal  length  as  shown  by  Hannan 
[24].  Consider  the  Cassegrain  antenna  in  figure  21.  The  main  reflector  is  a 
paraboloid  of  diameter  D and  focal  length  f,  and  the  subreflector  is  a hyper- 
boloid. The  real  focal  point  is  shown  between  the  reflectors,  and  the  equiva- 
lent single  paraboloid  (dashed)  has  a focal  length  f'  and  a diameter  D.  The 
equivalent  focal  length  can  be  written  as  [24] 


tan  (0/2) 

fl  = f ta-nTe;/^'* 


(55) 


For  a Gregorian  antenna,  the  subreflector  is  a concave  ellipse  as  shown 
in  figure  22.  The  equivalent  paraboloid  [24]  again  has  a diameter  D and  a 
focal  length  f'  given  by  eq  (55). 
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Dual  reflector  antennas  have  been  analyzed  by  using  physical  optics  on 
both  the  main  reflector  and  the  subreflector  [25].  It  was  found  that  the 
equivalent  paraboloid  approach  worked  well  for  both  Cassegrain  and  Gregorian 
antennas. 

4.5  Coax-to-Wavegu ide  Adapter 

In  this  section  we  analyze  the  out-of-band  response  of  a typical  coax-to- 
waveguide  adapter.  A knowledge  of  this  response  is  helpful  in  evaluating  the 
comparison  between  the  reflector  antenna  theory  and  the  experimental  results 
in  Section  3.3.  Such  adapters  are  also  used  in  many  other  microwave  antennas 
which  use  waveguide  feeds,  and  the  adapter  response  is  an  important  part  of 
the  total  antenna  system  response. 

The  probe  type  of  adapter  which  we  consider  is  shown  in  figure  23. 

Collin  [26]  has  analyzed  this  structure  using  a variational  technique,  and  he 
has  presented  some  numerical  results  for  in-band  frequencies  where  the  wave- 
guide supports  only  a single  propagating  mode.  Here  we  follow  Collin's  nota- 
tion and  formulation,  but  we  avoid  some  of  his  approximations  which  were  in- 
tended only  for  in-band  frequencies. 

The  main  task  is  to  compute  the  input  impedance  Z^n  of  the  probe  as  seen 
by  the  coaxial  cable  at  the  junction  (y  = 0).  Collin  [26]  has  derived  the 
following  double  integral  expression  for  Z,-n: 


where  S is  the  probe  surface,  _J_  is  the  surface  current  on  the  probe,  Q is  a 
dyadic  Green's  function  for  the  waveguide,  and  1^.  is  the  total  probe  current 
at  the  junction.  Collin  has  shown  that  eq  (56)  is  a variational  expression 
for  Z.jn;  thus  a first-order  approximation  to  J yields  a second-order 
approximation  to  Z^.  Consequently,  a sinusoidal  approximation  to  the  probe 
current  is  adequate,  and  J can  be  written  [26] 


Z. 


in 


jz  11  11  J (r)  * G ( r | r 1 ) • J (r ' ) da  da', 

t S S 


(56) 


J = a J sink  (d  - y), 
-y  o 


(57) 
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where  a is  a y-directed  unit  vector  and  JQ  is  an  arbitrary  constant.  We  as- 
sume that  the  current  distribution  is  uniform  in  the  circumferential  direction 
and  that  It  is  thus  given  by 


I.  = if  t J sin  kd, 
t o 


(58) 


where  t is  the  probe  diameter. 

In  evaluating  eq  (56),  it  is  convenient  to  write  Zin  as  a mode  sum  be- 
cause G is  a mode  sum. 
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where  Z 


nm 


Rnm  + ^nnr 


(59) 


Rnm  is  the  resistance  associated  with  the  nm  waveguide  mode,  and  Xnm  is  the 
reactance  associated  with  the  nm  waveguide  mode.  The  probe  excites  only  TM^m 
modes  that  are  transverse  magnetic  to  y.  The  fields  of  the  TM^  mode  have  n 
half  cycles  in  the  x direction  and  m half  cycles  in  the  y direction.  The 
terms  for  n even  do  not  contribute  to  Zn-n  because  the  probe  is  located  at  the 
center  of  the  waveguide  (x  = a/2).  For  propagating  modes,  the  propagation 
constant  3nm  is  given  by 

3 = j / k2  - ( n n/a ) 2 - (mn/b)2,  (60) 

nm 

where  a and  b are  the  waveguide  dimensions  indicated  in  figure  23.  For  evan- 
escent modes,  the  attenuation  constant  rnm  is  given  by 

r = / (mr/a)2  + (rWb)2  - k2.  (61) 

nm 

Only  the  propagating  modes  contribute  to  the  resistive  (real)  part  of  Zin,  and 
Collin  treated  the  in-band  case  where  only  the  TM^q  mode  was  propagating.  If 
we  apply  his  method  of  evaluation  to  an  arbitrary  nm  propagating  mode,  we  ob- 
tain the  following: 
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(62) 
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The  result  in  eq  (62)  agrees  with  Collin's  result  for  the  special  case,  n 
and  m = 0. 


= 1 


The  evanescent  modes  contribute  only  to  the  reactance  X^n.  For  suffi- 
ciently large  values  of  m,  we  follow  the  Poisson  summation  formula  of  Collin 
and  obtain  the  following: 
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n=l,3,< 


nm 
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K 0<  t/2), 

o m 


(63) 


where  k = / (mrr/b)2  - k2  , 
m 

and  KQ  is  a zero-order,  modified  Bessel  function  [14].  The  result  in  eq  (63) 

is  valid  for  m > M where  k^  a >>  1.  The  sum  on  m from  M to  00  converges  rapid- 

ly because  KQ  decays  exponentially  for  large  argument  [14].  For  m < M,  we 
perform  the  n summation  using  the  same  asymptotic  technique  that  Collin  used 
for  m = 0.  The  technique  is  based  on  the  following  asymptotic  behavior 
of  r for  large  n: 

Jim  r -1  ~ -i  -V2  ( -±  )2  [(  STL  )-  k2l  . (64) 

n+°°  nm  nu  2 nn  ' L ^ b > J 


However,  we  must  include  a term,  exp  (-  rnm^)»  which  Collin  dropped  because  it 
was  exponentially  small  for  in-band  frequencies.  When  this  is  done,  the  sum 
for  m < M is  given  by 
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The  primes  on  the  summations  indicate  the  omission  of  the  propagating  modes 
that  are  computed  by  eq  (62)  and  the  P's  on  the  final  sum  indicate  only  the 
propagating  modes.  The  terms  in  the  infinite  double  sum  go  to  zero  as  n“J  for 
large  n because  of  the  asymptotic  behavior  of  rnm~'*'  as  given  by  eq  (63). 

A computer  code  for  n has  been  written  and  the  double  summation  over  n 
and  m is  performed  as  indicated  in  figure  24.  The  boundaries  between  the 
three  regions  depend  on  the  parameters  a,  b,  and  k,  but  the  computer  code 
determines  the  appropriate  boundaries  automatically.  The  computer  code  is 
fast  because  both  the  n summation  in  eq  (65)  and  the  M summation  in  eq  (63) 
converge  rapidly.  The  code  has  been  checked  against  Collin's  curves  for  X^n 
and  R.jn,  and  agreement  has  been  obtained  to  graphical  accuracy. 


We  now  consider  the  transmission  of  power  from  the  coaxial  cable  to  the 
waveguide.  We  assume  that  the  coaxial  cable  supports  only  a dominant  TEM  mode 
and  that  the  coaxial  cable  has  a real  characteristic  impedance  Rc.  Then  the 
voltage  reflection  coefficient  ry  is  given  by 


r 


V 


(66) 


The  power  reflection  coefficient  is  then  given  by  | rv | 2 , and  the  transmission 
coefficient  for  the  total  power  supplied  to  the  waveguide  T^  is 
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(67) 


We  can  also  define  the  transmission  coefficient  for  the  nm  mode  Tpm  as  the 
ratio  of  the  power  transmitted  to  the  nm  mode  to  the  incident  power  in  the 
coaxial  cable.  By  using  the  relationship  between  base  current  It  and  Rnm 
[26],  Tpm  is  found  to  be 


4 R R 
nm  c 

nm  (R.  + R )2  + X.  2 * 

in  c in 


(68) 


From  eqs  (67)  and  (68),  it  is  easy  to  see  that  the  total  power  transmitted  to 
the  waveguide  is  equal  to  the  sum  of  the  power  transmitted  to  the  individual 
propagating  modes: 
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(69) 
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In  the  receiving  antenna  problem,  we  are  interested  in  the  reciprocal  problem 
where  a propagating  waveguide  mode  transmits  power  to  the  coaxial  cable.  By 
reciprocity,  it  can  be  shown  that  the  power  transmission  coefficient  Tnm  as 
given  by  eq  (68)  also  applies  to  this  case. 

A computer  code  was  written  to  compute  Tnm  from  eq  (68),  and  it  was  ap- 
plied to  an  S-band  coax-to-waveguide  adapter.  The  adapter  parameters  were 
chosen  to  match  those  used  in  the  reflector  antenna  which  was  studied  experi- 
mentally by  Cown,  et  al.  [5,27]:  a = 7.112  cm,  b = 3.302  cm,  d = 1.9  cm,  and 

i = 2.4  cm.  The  radius  of  the  probe  was  not  given,  but  we  assumed  a value  of 
t/2  = 3.5  mm  in  order  to  yield  a small  reactance  at  3 GHz.  The  assumed  value 
of  t is  not  critical  because  the  results  are  only  weakly  dependent  on  the 
probe  radius.  The  transmission  coefficients  as  a function  of  frequency  are 
shown  in  figures  25  and  26.  Mote  that  the  transmission  coefficient  for  the 
dominant  TM^q  mode  is  nearly  unity  from  about  2.5  to  5.0  GHz.  Above  5 GHz, 
higher  order  modes  begin  to  appear,  and  the  results  are  very  frequency 
sensitive.  Some  related  calculations  have  been  done  by  Cown  and  Ryan  [27], 
but  they  did  not  compute  the  reactance  of  the  probe  which  is  required  to  com- 
pute transmission  coefficients  as  indicated  by  eq  (68).  At  10  GHz,  there  are 
five  propagating  modes  (TM^q,  TM^,  TM^q,  TM^,  and  TM^)  which  are  coupled  to 
the  adapter.  In  addition,  there  are  propagating  TE^m  (transverse  electric) 
and  propagating  TM^m  modes  with  n even  [27],  but  these  modes  do  not  couple  to 
the  adapter.  All  of  the  transmission  coefficients  have  a null  at  7.9  GHz 
because  the  sinusoidal  current  assumption  in  eq  (57)  yields  an  infinite  value 
of  Zn-n  at  kd  = tt. 

In  the  earlier  comparison  between  theory  and  experiment  for  on-axis  gain 
in  figure  14,  the  experimental  values  were  well  below  the  theoretical  re- 
sults. The  main  problem  with  the  comparison  was  that  the  theory  included  only 
the  antenna  response  while  the  experiment  included  the  antenna  plus  the 
adapter  response  [27],  No  attempt  was  made  to  separate  the  antenna  and 
adapter  responses  because  the  out-of-band  response  of  the  adapter  was  not 
known.  Even  though  we  now  have  theoretical  results  for  the  adapter  response, 
it  is  not  clear  how  to  combine  the  results  with  the  reflector  antenna  theory 
at  out-of-band  frequencies  where  higher  order  modes  are  propagating  in  the 
waveguide.  If  we  take  the  total  transmission  coefficient  Tj.  as  given  by  eqs 
(67)  or  (69)  and  multiply  by  the  effective  aperture  of  the  antenna,  then  a 
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corrected  frequency  response  is  obtained  as  shown  in  figure  27.  The  correc- 
tion can  be  made  in  other  ways,  but  no  precise  correction  can  be  made  because 
the  antenna  theory  gives  only  the  total  waveguide  power,  not  the  modal  con- 
tent. However,  the  adapter  correction  does  improve  the  agreement  significant- 
ly, particularly  below  7.5  GHz.  The  null  in  the  adapter  response  at  7.9  GHz 
is  not  apparent  in  the  experimental  data,  but  no  measurements  were  made  be- 
tween 7.5  and  8.0  GHz.  Above  8 GHz,  the  agreement  is  not  as  good,  but  this  is 
probably  because  the  adapter  theory  is  less  accurate  when  the  probe  is  longer 
than  a half  wavelength. 

It  might  be  possible  to  improve  the  adapter  theory  for  higher  frequencies 
(kd  > it)  by  using  a more  accurate  distribution  for  the  probe  current  [26],  but 
it  is  doubtful  that  the  detailed  structure  above  8 GHz  in  figure  27  could  be 
obtained.  The  thin  wire  and  small  coaxial  gap  assumptions  are  not  really 
valid  above  8 GHz,  and  also  the  constructive  and  destructive  interference  of 
the  higher  order  waveguide  modes  is  probably  responsible  for  the  rapid  fre- 
quency variations  in  figure  27.  The  main  point  of  figure  27  is  that  the 
theory  matches  the  experimental  results  fairly  well  in  the  region  where  the 
adapter  theory  is  valid  (below  7.5  GHz),  and  we  can  expect  the  reflector 
antenna  theory  to  remain  valid  for  higher  frequencies  because  it  is  a high 
frequency  theory.  It  might  be  desirable  to  study  the  isolated  adapter  re- 
sponse experimentally,  but  it  might  be  difficult  to  perform  a meaningful  ex- 
periment because  of  the  higher  order  modes. 


5.  Conclusions 

The  response  of  reflector  antennas  to  out-of-band  (above-band)  frequen- 
cies has  been  studied  in  a two-step  analysis.  In  Section  2,  the  electric  and 
magnetic  fields  and  the  Poynting  vector  in  the  focal  region  of  a paraboloidal 
reflector  have  been  determined  by  a physical  optics  integration.  The  large 
f/D  approximation  is  found  to  be  sufficiently  accurate  for  most  realistic 
antennas  (f/D  > 0.3),  and  numerical  integration  is  therefore  not  required. 

The  second  step  of  the  analysis  is  an  integration  of  the  Poynting  vector  over 
the  aperture  of  the  feed  horn  to  obtain  the  received  power.  As  shown  in 
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Section  3,  this  integration  can  also  be  done  by  an  analytical  approximation  in 
order  to  eliminate  the  need  for  numerical  integration.  The  results  yield  a 
generalized  effective  aperture  (the  total  received  power  divided  by  the  inci- 
dent power  density)  for  the  antenna,  and  a fairly  simple  expression  gives  both 
the  frequency  response  and  the  receiving  pattern. 

The  theoretical  results  have  been  compared  with  published  measured  re- 
sults for  a symmetrical  paraboloid  with  f/D  = 0.32  over  a frequency  range  from 
3 GHz  (in-band)  to  10  GHz.  The  theoretical  pattern  shape  matches  the  envelope 
of  the  measured  pattern,  but  the  measured  gain  falls  well  below  the  theoreti- 
cal result.  When  the  transmission  characteristics  of  the  coax-to-waveguide 
adapter  are  taken  into  account,  the  agreement  between  theory  and  experiment  is 
much  better.  Two  important  results  of  the  theory  are  that  the  beamwidth  at 
out-of-band  frequencies  is  approximately  the  same  as  the  in-band  beamwidth  and 
that  there  are  no  high  sidelobes.  These  results  are  in  agreement  with  the 
published  measured  data  [27]. 

In  Section  4,  a number  of  extensions  to  the  basic  physical  optics  model 
are  discussed.  The  level  of  the  distant  sidelobes  is  examined  using  GTD  and 
is  found  to  be  in  approximate  agreement  with  physical  optics.  The  transient 
fields  in  the  focal  region  are  examined  for  pulse  excitation,  and  the  pulse 
stretching  is  in  agreement  with  simple  geometrical  arguments.  The  effect  of 
reflector  roughness  is  examined,  and  the  reduction  in  gain  or  effective  aper- 
ture can  be  calculated  if  the  surface  tolerance  is  known.  Extensions  of  the 
theory  for  a symmetrical  reflector  to  the  cases  of  offset  or  dual  reflectors 
is  accomplished  by  using  equivalent  focal  length  and  diameter.  Also,  the  pre- 
vious theory  for  the  response  of  a coax-to-waveguide  adapter  is  extended  from 
the  in-band  case  [26]  to  the  out-of-band  case.  This  theory  might  be  useful 
for  other  microwave  antennas  using  the  same  type  of  adapter  in  the  feed  sys- 
tem. 

Other  extensions  and  improvements  to  the  model,  such  as  a more  precise 
treatment  of  the  feed  horn,  are  certainly  possible,  but  the  physical  optics 
model  presented  here  is  probably  adequate  for  typical  out-of-band  applica- 
tions. The  analysis  of  other  common  microwave  antenna  types,  such  as  phased 
arrays,  is  recommended  as  having  highest  priority. 
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Appendix  A 

Evaluation  of  the  <t>i  Integration 


The  first  step  in  the  evaluation  of  eq  (6)  is  to  write  the  phase  term  in 
the  following  form: 


ns  + R = Mo  ~ M1  cos^l" 


(Al) 


where  MQ,  M^,  and  ij>  are  independent  of  4>]_.  From  eqs  (7)  and  (8),  we  find  that 
the  three  quantities  are  given  by 


M = A + 
o o 


5 


M = / (A  + B )2 * *  + B 2 

1 CCS 


(A2) 


and  = sin~^  » 

where  AQ  = (f  - rj)  coses. 


Ac  = rj.  sinO^  sin0$. 


Bc  = - r2  sin0j  cos^^ 
and  = - r2  sin0^  sin<t>2. 

The  next  step  is  to  express  the  amplitude  terms  in  eq  (6)  in  terms  of  a 
Fourier  series  in  <)>]..  When  this  is  done,  it  is  found  that  the  series  trun- 
cates at  n = 2,  and  dE_  and  d]l_  can  be  written: 

dE  = - (Mo  + M1  cos4,). 

I \ [Cxncos  n+l  + Sxns1n  n+l>  u-x  + \ [V°S  "*1  + Syns1n  n*l]  h 

n=0  n=0  J j j 

2 2 sin0, 

+ Zn  tCzncos  "*1  + Szns1n  n *1>  -uz>  1 + cose.  d0> 

n=0  1 
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and 


(A3) 


E.  f 


dH 


_•  _i -jk  (M  + M,  cos^) 

" J > c U 1 • 


1 [Dxncos  n*l  + Fxns1n  n*l]  -x  + K CV°S  "*1  + Fyns1n  n*l]  u- 
n=0  n=0  J J 


2 sine. 


+ s.  [DznC0S  n*l  + Fznsin  n4,l]  a2l  l + cos'e,  de’  d<fl 
n=U  I 


The  Fourier  coefficients  are  obtained  from  the  terms  inside  the  brackets  in  eq 
(6)  and  found  to  be: 


sin2e, 

C ~ = a (cosQ-i  + V?  t — t r— ).  S n = 0, 

xO  y v 1 2 1 + cos0^;*  xO  * 


Cxl  ■ 0 
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Cx2  " 1/2 
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sin20. 
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From  Fourier-Bessel  expansions  [2.7],  we  can  derive  the  following  useful 
integral : 


I | cos  n 4>  . -jk  M,  cos(  <*», 

1 1 sin  nct>  ; 


^^d<t>!  = 2tt  jn 


{ cos  nif  } J (k  M,), 
1 sin  m|i  J n v V * 


(A5) 


where  Jn  is  the  nth  order  Bessel  function.  By  using  eqs  (A5)  in  (A3),  the 
integrations  can  be  carried  out  to  yield  the  desired  expressions  for  E_  and  H_ 
in  eg  (13). 
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Appendix  B 

Open-ended  Waveguide 


Consider  the  two-dimensional,  open-ended  waveguide  of  width  w in  figure 
28.  A plane  wave  is  incident  in  the  xz  plane,  and  the  incident  electric  field 
is  y polarized: 


E , = E e 
yi  o 


jk  (z  cose  + x sin9) 


(Bl) 


where  0 is  the  incidence  angle  measured  from  the  z axis, 
the  incident  Poynting  vector  Sz  is  given  by 


The  y component  of 


E I2 
o| 

2n 


cos  0 . 


(B2) 


From  eg  (23),  the  approximate  received  power  per  unit  length  P is  given  by 


P = Pn  COS0, 
r o * 


( B3) 


where  PQ  = w SQ  and  S 


2n  ‘ 


If  we  make  the  Kirchhoff  approximation  in  the  waveguide  aperture  (z  = 0, 
x|  < w/2),  then  the  electric  field  is  given  by 


'y  U_ 


= E e 


jkx  sin9 


( B4 ) 


z=0 


A general  waveguide  mode  expansion  for  the  electric  field  E^  is 


r _ v r (2  n-1)  ttx  _J'kxez 

Ey  - I lcn  cos  e 


2n  ttx  jkxoZ! 

+ sn  S1n  w e f 


(B5) 


■ . . . o r (2  n"^)  77  i2 

where  k = / k2  - 

xe  w 

and  k = / k2  - ( ) . 

xo  w 

The  cutoff  wavelength  Ac  is  determined  by  setting  kxe  = 0 for  the  even  (in  x) 
modes  and  by  setting  kxo  = 0 for  the  odd  (in  x)  modes: 
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Even: 


= 0 
xe 

, 2w 

c 2 n-1  » 

X 

o 

II 

o 

X = - . 
c n 

( B6) 


The  unknown  coefficients  are  determined  by  matching  the  electric  field  at  y = 
0 in  eqs  (B4)  and  (B5) 


E ejkx  s1ne  = l [c  cos  ■n-1)  ” + s 

L n ui  n ui  -* 


n=l 


(B7 ) 


By  using  the  orthogonal ity  relationships  for  the  cosine  and  sine  functions  in 
eg  (B7),  we  obtain  the  following  expressions  for  cn  and  sn: 


c = En 
n o 


r kw  . „ . (2  n-1)  tt 
sin  [ sine  + A — ] 

^ r kw  . ~ , (2  n-1)  tt 
[ — 2"  sine  + — 


] 


and 


sin  [-^  sine  - -y-— -] 

sine  - 


(B8) 


s_  = 


r kw  sine  , ni^ 

sin  [ 2 + ~-i 

r kw  sine  nn-. 

L 2 2J 


. rkw  sine  nir-, 

sin  [ 2 2-* 

rkw  sine  nn-. 

L 2 " 2J 


The  sin  ( )/(  ) in  eq  (B8)  have  peaks  where  the  transverse  wave  numbers  of  the 
even  and  odd  waveguide  modes,  (2  n-ljn/w  and  2mr/w,  match  the  transverse  wave 
number  of  the  incident  field,  k sine. 

In  order  to  compute  the  power  in  the  waveguide  modes,  we  need  an  expres- 
sion for  the  magnetic  field  H: 


H = — V x u.  E... 


GJU 


-y  y 


From  eq  (B9)  the  magnetic  field  components  Hx  and  Hz  are 


1 


3 E 1 d L 

H = -r^ and  H = -? 

x jwp  3 z z 1‘-"‘ 


z 
3E 

jwp  3x* 


(B9) 


(BIO) 


The  power  per  unit  length  transmitted  by  the  waveguide  Pg  is  given  by 

w/2 

P = / Re  (E  H *)  dx. 

9 -w/2  y z 


(Bll) 
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By  substituting  eqs  (B5)  and  (BIO)  into  (Bll)  and  carrying  out  the  x integra- 
tion, we  obtain 


P 

9 

where  Pne  = 4F7  lcrj2  Re  (kxe> 
and  Pno  = 4R  lsnl2  Re  (kxo>- 


00 


z 

n=l 


+ 


P 

no 


), 


(B12) 


For  modes  above  cutoff,  kxe  and  kXQ  are  real.  For  modes  below  cutoff,  kxe  and 
kx0  are  purely  imaginary.  Thus,  only  the  propagating  modes  carry  power,  and 
the  summation  in  eq  (B12)  can  be  truncated. 


A check  on  the  accuracy  of  the  Kirchhoff  approximation  for  the  aperture 
field  in  eq  (B4)  can  be  made  by  comparing  the  waveguide  power  Pg  in  eq  (B12) 
with  the  received  power  approximation  Pp  in  eq  (B3).  In  general,  the  evalua- 
tion of  eq  (B12)  must  be  done  numerically  for  specific  parameters.  However, 
the  limit  for  large  kw  can  be  evaluated  analytically  for  on-axis  incidence, 
0=0.  In  this  case,  the  mode  coefficients  are  given  by 


= 4 (-D 


n-1 


cn  (2  n-1)  tt  E( 


and  s„  = 0. 


(B13) 


For  large  kw,  kxe  « k,  and  Pn  is  given  by 


P = 
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4w 

E 

0 

2 

TT" 

- n 

Z (2  n-1)'2. 
n=l 


The  infinite  sum  has  a known  result  [28]: 


l (2  n-1) 
n=l 


-2 

8* 


Substituting  eq  (B15)  into  eq  (B14),  we  have 


w|E 


P = 
9 


2n 


■ W So  * Pr* 


(B14) 


(B15) 


(816) 


Thus  Kirchhoff  approximation  yields  a consistent  result  for  the  received  power 
and  the  waveguide  power  when  kw  is  very  large. 
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In  figures  29-31,  we  show  results  for  the  received  power  and  the  wave- 
guide power  as  a function  of  0 for  various  values  of  w/x.  In  all  cases,  the 
various  powers  are  normalized  to  the  incident  power  PQ  at  0=0.  As  w/X  is  in- 
creased from  1.5  to  10,  the  agreement  between  Pr  and  Pg  improves  for  all  inci- 
dence angles.  In  all  cases,  the  powers  in  the  first  even  and  odd  modes  are 
shown  for  comparison.  As  w/x  is  increased,  the  mode  patterns  show  more  rapid 
oscillation,  and  the  number  of  propagating  modes  increases.  For  w/X  = 1.5, 
there  are  only  2 propagating  modes,  but  for  w/X  = 10,  there  are  19  propagating 
modes. 
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Figure  1.  Geometry  for  a plane  wave  incident  on  a symmetrical  parabolic 
ref 1 ector. 
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Figure  2.  Normalized  Poynting  vector  for  on-axis  incidence. 
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Figure  3.  Normalized  Poynting  vector  for  off-axis  incidence. 
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Figure  4.  Geometry  for  a circular  feed  horn. 
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Figure  5.  Normalized  Poynting  vector  and  received  power  for  on-axis 
incidence. 
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Figure  6.  Approximate  integration  surface  (dashed)  for  off-axis  incidence. 
The  geometrical  optics  point  is  located  inside  the  circle. 
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Figure  7.  Approximate  integration  surface  (dashed)  for  off-axis  incidence. 
The  geometrical  optics  point  is  located  outside  the  circle. 
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Figure  8.  Normalized  effective  aperture  for  off-axis  incidence. 
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Figure  9.  Comparison  of  analytic  and  numerical  results  for  the  effective 
aperture. 
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Figure  10.  Comparison  of  theory  and  experiment  [3.3]  for  the  in-band 
frequency  of  3 GHz. 
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Figure  11.  Comparison  of  theory  and  experiment  [5]  for  an  out-of-band 
frequency  of  6 GHz. 
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Figure  12.  Comparison  of  theory  and  experiment  [5]  for  an  out-of-band 
frequency  of  8 GHz. 
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Figure  13.  Comparison  of  theory  and  experiment  [5]  for  an  out-of-band 
frequency  of  10  GHz. 
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Figure  14.  Comparison  of  theory  and  experiment  [27]  for  gain  as  a function  of 
f requency. 
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Figure  15.  Geometry  for  diffraction  from  the  edge  of  the  reflector. 
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Figure  16.  Geometry  for  direct  illumination  of  the  feed  horn. 


56 


0 


-10  - 


-20  - 


D-1.22  m 
P - 3.39  cm 
f/D  = 0.32 
freq.  = 6 GHz 


Phys.  Opt. 


Direct  Ilium. 


1 20 


Figure  17.  Various  contributions  to  the  wide  angle  receiving  pattern. 
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Figure  18.  Normalized  impulse  response. 
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Figure  19.  Roughness  loss  as  a function  of  normalized  frequency. 
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Figure  20. 


Geometry  for  an  offset  paraboloid. 
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Figure  21.  Geometry  for 
(dashed)  has  a focal 


a Cassegrain 
length  f'. 


antenna. 


The  equivalent  paraboloid 
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Figure  22.  Geometry  for  a Gregorian  antenna.  The  equivalent  paraboloid 
(dashed)  has  a focal  length  f'. 
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Figure  23.  Probe  type  of  coax-to-wavegu ide  adapter. 


63 


Figure  24.  Computational  scheme  for  evaluating  Z1n.  The  nm  modes  are  grouped 
into  three  regions:  propagating  modes,  evanescent  modes  for  m < M,  and 
evanescent  modes  for  m > M. 
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Figure  25.  Transmission  coefficients  for  an  S-band  coax-to-wavegu ide 

adapter.  The  cutoff  frequencies  for  the  dominant  ( 0 ) and  first  higher- 
order  (TM^)  modes  are  indicated  by  arrows. 
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Figure  26.  Transmission  coefficients  for  the  higher  order  modes  of  an  S-band 
coax-to-waveguide  adapter.  The  cutoff  frequencies  for  the  higher  order 
waveguide  modes  are  indicated  by  arrows. 
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Figure  27.  Comparison  of  theory  (with  adapter  correction)  and  experiment  [27] 
for  antenna  gain  as  a function  of  frequency. 
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Figure  28.  Plane  wave  incident  on  an  open-ended,  parallel-plate  waveguide. 
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Figure  29.  Total  received  power  and  the  contributions  from  the  n=l  even  and 
odd  modes.  Only  two  modes  are  propagating. 
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Figure  30.  Total  received  power  and  the  contributions  from  the  n=l  even  and 
odd  modes.  Five  modes  are  propagating. 
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Figure  31.  Total  received  power  and  the  contributions  from  the  n=l  even  and 
odd  modes.  Nineteen  modes  are  propagating. 
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